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FOREWORD 


While  tracking  a  target  it  is  desirable  to  know  how  well  the  target 
tracking  algorithm  thinks  it  knows  the  target’s  location.  A  common  method  used 
to  indicate  the  quality  of  the  track  is  to  display  the  one-sigma  uncertainty 
ellipse.  This  indicator  of  the  track’s  quality  is  derived  directly  from  the 
tracker’s  covariance  matrix. 

The  following  are  two  methods  for  deriving  the  one-sigma  uncertainty 
ellipse.  The  first  method  uses  the  trigometric  equations  of  an  ellipse,  and  the 
second  uses  matrix  notation  to  accomplish  the  same  task.  The  last  part  of  this 
paper  describes  a  parametric  method  of  plotting  any  ellipse. 


Approved  by: 


C.  A.  KALIVRETENOS ,  Head 
Sensors  and  Electronics  Division 
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CHAPTER  1 

DERIVATION  OF  THE  ONE-SIGMA  PROBABILITY  ELLIPTIC  CONTOURS 
USING  TRIGONOMETRIC  EQUATIONS 


GIVEN:  A  Gaussian  joint  probability  density  function: 


f(x,y) 


1 


2(1-P2) 


2 

2pxy  +  y_ 

a  a  a 

x  y  y 


The  one-sigma  contour  is  defined  by 


2Cl-p2) 


a  a 

x  y 


=  1. 


(1) 


We  want  to  find  the  semi-major  axis,  the  semi-minor  axis,  and  the  angle  of 
inclination  of  the  one-sigma  contour  equation. 

The  angle  of  inclination,  0,  will  first  be  found.  To  do  this  we  rotate  the 
axis  of  the  ellipse  until  the  xy  coefficient  vanishes.  Rewriting  Equation  (1)  as 


2  2  2  2  2  2  2 
a  x  -  2a  xy  +  a  y  =  2(l-p  )  a  a  ,  where  a  =  pa  a  , 
y  xy  x  x  y  ’  xy  x  y 


and  substituting  in  the  equations  of  rotation,  from  Figure  1,  i.e., 


x  =  x*cos0  -  y*sin0 


y  =  x*sin0  +  y*cos0, 
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\ 


\ 


\ 


X*  cos© 

FIGURE  I.  TWO-DIMENSIONAL  ROTATION  THROUGH  AN  ARBITRARY  ANGLE,  0 


results  in 


2  2  2  2  *2 
(a  cos  9  -  2a  sin9cos0  +  a  sin  0)x 
y  xy  x 

2  2  2  2  * 

+  (2a  sin  0  +  2(a  -  a  )sin0cos0  -  2a  cos  0)x  y 

xy  x  y  xy  J 

2  2  o  2  *  2 

+  (cr  cos  6  4*  2a  sinBcosB  +  a  sin  B)y  =  2(1  -  p 
x  xy  y 


To  determine  0  we  need  only  to  set  the  x  y  coefficient  to 
the  angle  of  rotation  needed  to  make  the  major  axis  of  the 
x  -axiSo  If  the  identities 


jo  o 

x  y 

zero;  i.e.,  0  is 
ellipse  fall  on  the 


cos^B  -  sin^B  **  cos20 
2sinBcosB  -  sin2B 


(2) 
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are  substituted  into  the  x  y  coefficients  of  Equation  (2)  and  then  this 
coefficient  is  equated  to  zero,  the  result  is: 


2  2 

(a  -  a  )  sin29  -  2 a  cos29  =  0. 
x  y  xy 


(3) 


Thus  the  angle  we  must  rotate  the  axes  is 


2o 

0  =  1/2  ARCTAN  — — \ 

o  -a 

x  y 


To  find  the  semi-major  axis  and  semi-minor  axis  we  substitute  the 
identities, 


2n  1  +  cos29 
cos  9  - 


2„  _  1  +  cos29 
sin  9  =  - 


2sin0cos9  =  sin29 , 


into  Equation  (2)  resulting  in 


2  /I  +  cos29 \ 


.  .  2  (l  -  cos2e\  *2 

o  I - x - 1  -  a  sm2e  +  a  - ~ - lx 

y\2/xy  x  \  2  / 


2  2 

a  (1  -  cos20 )  +  (a  -  a  )  sin20  -  a  (1  +  cos2@) 
xy  x  y  xy 


2/1  +  cos20\  ,  .  „„  ,  2/1  -  cos2e\ 

( - 2 - )  +  a _  sm20  +  ^ -  ) 


xy 


2  / 


*2 


*  * 
x  y 


(4) 


2,2  2 


y  -  2(l-p  )ax  ay 


If  we  now  write  Equation  (3)  as 


2a 


tan  20  - 


xy 


2  2 

a  -  a 
x  y 
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and  depict  this  equation  in  Figure  2,  we  can  see  the  relationships 


2a 


sin29  = 


xy 


2  2 

a  -  a 


cos 


29  - L 


where 


Y  =  \(a  2  -  a  2)  +  4a  2 

x  y  xy 


From  these  relationships  and  Equation  (4)  the  x  coefficient.  A,  becomes 
A  =  1/2 

|_  y  x  y  x 

2 


2  2  2  2 

a  +  a  +  (a  -  a  )  cos29  -  2a  sin20 
y  x  y  x  xy 


1/2 


=  1/2 


2  2 

9  9  |°  -  a 

a  2  +  a  2  ^  x 

l  y  x 


4  a  “ 

Ll _ EL 

Y  Y 


2  2 

a  +  a 

y  x 


■i 


2  2 1  .  2 
-  a  I  +  4a 
x  y  /  xy 


Likewise  the  y  coefficient,  B,  is 
B  =  1/2 


2  2 
a  +  a  + 

y  x 


2,2  ,  ,  2 

a  )  +  4a 

y  xy 


Also  we  can  rewrite  the  term  under  the  radical  as: 

,  2  2x2  .  2 

Ca  -  a  )  +  4a 

x  y  xy 

4  _  2  2  4  ,  2  2  .2  2  .  2 

=  a  -2a  a  +a  +4aa  -4a  a  +4a 
x  xyy  xy  xy  xy 


-  (a  +  a  )  -  4a  ^a  (1  -  p  ),  where  p  =  — . 

x  y  xy  aa 

x  y 
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FIGURE  2.  GEOMETRIC  INTERPRETATION  OF  EQUATION  3 


+  4  a 


2 


xy 
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2  2 

Letting  b  =*  a  +  a  we  find 
x  y 


A  =  1/2  jb  -  |V  -  4ox2ay2(l  -  p  2) 


i"’} 


and 


B  =  1/2 


{b+  [b2-4<,_Vu  -P2) 


x  y 


»]■'*[. 


Now  we  can  rewrite  Equation  (1),  rotated  9  radians,  as 


*2  *2  2  2  2 
AxZ  +  ByZ  =  2(1  -p  ;  a  a  Z 

x  y 


This  equation  can  be  equated  to  the  standard  elliptic  equation  as  follows: 


o  f  i  2  v  2  2 

2(1  -  p  )ax  ay 


*2  , 
x  + 


2  2 

_ B  *2  _  .  x  .  / 

2  2  2  ^  ^  2+2 
2(1  -p  )axayZ  a2  c 


From  this  we  see  the  semi-major  axis  is 


a  = 


v 


,,,  2,  2  2 
4(1  -  p  )  a  a 
x  y 


Frw^n, 


and  the  semi-minor  axis  is 


c  = 


w  -  2  v  2  2 

40  ~  p  )  a  x  ay 


H‘j-sVd-!>n 
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CHAPTER  2 

DERIVATION  OF  THE  ONE-SIGMA  PROBABILITY  ELLIPTIC  CONTOUR  USING  MATRICES 


GIVEN:  A  Gaussian  joint  probability  density  function 


f(x,y)  = 


2tt  |  Q  | 


1/2 


-1/2(XTQ  1X) 


where  X  = 

and  the  matrix  Q  ^  is  the  inverse  of  the  covariance  matrix,  i.e., 


x-1 


2  2f 1  -  2^ 

av  a„  U  p  ) 


x  y 


-pa  a 

x  y 


-pa  a 

x  y 


The  one-sigma  contour  is  defined  by 
1/2XTQ  “*X  -  1. 

Since  Q  is  a  Hermitian  matrix  of  rank  2  a  substitution  can  be  made 
Q-1  =  PDPt 

where  D  is  a  diagonal  matrix 


D  = 


\  0 


0  K 


^9  A  ■ 


P  is  a  matrix  formed  from  the  eigenvectors  of  Q  ^  corresponding  to  the 
eigenvalues  ^  and  i.e., 


P  = 


cos0  -sin0' 
sin0  cos0 


(5) 

(6) 

(7) 
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0  represents  an  angle  of  rotation  between  the  x,y  coordinate  system  and  another 
system  oriented  in  the  directions  of  ellipse’s  major  and  minor  axes.  We 
calculate  the  eigenvalues  as  follows: 


|Q  1  -  All 


=  det 


2  2  2 
a  a  ^(l-pz) 
x  y 


-  A 


-pa  a 

_ x  y 

2  2m  2^ 
a  a  U-P  ) 


-pa  a 

x  y 

a  2 a  2(1-P2) 

x  y 


2  1(^  1\ 
°x  CTy  (1~p  } 


=  0  . 


A 


Solving  for  the  eigenvalues  A \  and  A2  we  get 

2 


2  ,  2 

a  +  a 

a  =  _* _ l_ 

1 


-  ( 


,  2  .  2,  .12,.  2. 

(a  +  a  )  -4a  a  (1  -  p  ) 
x  y  x  y 


1/2 


22  2 

2ax  °y  (1  "  p  > 


X2  = 


[‘ 


2  2 . 


a  +  a  +  |(a  +  a  )  -  4a  a  (1  -p  ) 

”  Y _ x  y 


1/2 


-5  2  2/n  2n 

2CTX  °y  (1  "  P  ) 


To  show  the  relationship  between  the  eigenvalues  and  the  major  and  minor  axes  of 
the  ellipse,  the  coordinates  are  rotated  by  use  of  the  P  matrix.  From  Equation 
(6) 

XtQ-1X  =  2 

and,  using  Equation  (7)  X^PDP^X  =  2,  let 


X*  = 


where 


X  =  PX" 


Then  XT  =  X*TPT  and  PT  =  P-1 . 

Also,  it  follows  that  X*^DX*  =  2 

and  Ajx*2  +  A2y*2  =  2. 


Equating  this  equation  to  the  standard  equation  for  an  ellipse  we  get: 
*2  *2  *2  *2 
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So  the  semi-major  axis  is 


a  = 


,22. .  2, 

Oy  (l  -  P  ) 


-  [b2  -  40..V.2  (1  -  P2)l  1/2  ’ 


x  y 


and  the  semi-minor  axis  is 


c  - 


,  2  2  ..  2, 
4ax  °y  (1  “  p  ) 


2  \  b  + 


172  ,  2  2  ,  2, 

b  -  4o  _  a  w  (1  -  p  ) 


x  y 


T72  * 


The  angle  of  inclination  of  the  ellipse  is  obtained  by  solving  for  9.  Since 
Q-1  -  PDPt, 

then  D  =  P-^Q-lp.  Now  we  may  rewrite  this  equation  as: 


T  -1 
P  Q  XP  = 


„  2n  2  n  2 
a  a  -  a 

x  y  xy 


“cos0 

sin0~ 

_  L. 

a 

y 

-a  “1 

xy 

COS0 

-sin0 

-sin9 

COS0_ 

-a 

L  xy 

2 

ax  J 

sin0 

COS0 

r  2  2 


2  2  2 

a  o  -  a 

x  y  xy 


2  .  2 


2  2, 


a  cos  0  -  2a  sinecose  4-  a  sin  0  a  (sin  0-cos  0)  +  (a  -a  )sin0cos0 
y  xy  x  xy  x  y 

2„  2_ 


2  2, 


2  ,  2 


a  (sin^0-cos^0 )  4-  (a  ^)sin0cos0  a  ^sin^0  4-  2a  sin0cos0  +  a  cos  0 

~  y  xy  x  _ 


I-  xy  '  ' "  x  ~  y 

where  we  made  use  of  the  notational  change 

pa  a  =  a 
x  y  xy  • 

Setting  the  off  diagonal  terms  equal  to  zero  yields 

2CT  XV 

tan2e  =  -  2  _7  2 

ax  Oy 

and  9  =  j  ARCTAN  — — $ 

ax  ay 

The  angle  of  inclination,  9 ,  is  defined  as  the  counterclockwise  rotation  from 
the  x  axis  to  the  major  axis  of  the  ellipse. 
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CHAPTER  3 

PLOTTING  A  GENERAL  ELLIPSE  WITH  ITS  CENTER  LOCATED  AT  (h,  k) 


The  equations  for  a  general  ellipse  with  inclination,  0  ,  and  centered  at 
(h,  k)  can  be  written  as 


£x*  y*  lj  =  £a  cos  <j>  b  sin  <J>  lj 


~cos0  sin0  0 

-sin  0  cos  0  0 

h  k  1« 


or 


x 


* 


-  a  cos  <f>  cos  0  -  b  sin<()  sin  0  4-  h 


y  -  a  cos  (j)  sin  0  4*  b  sin  (j)  cos  0  4*  k  . 

Here  a  is  the  semi-major  axis  in  the  x  direction,  b  is  the  semi-minor  axis  in 
the  y  direction  and  $  is  a  variable  for  plotting  that  allows  us  to  plot  the 
elliptic  contour. 

If  we  let  the  angle,  ,  increment  in  a  A(j)  step,  the  next  location  for  the 
locus  in  the  x*y*  plane  will  be 


* 

x  =  a  cos  (  (j)  +  A<j> )  cos  0  -  b  sin  (<j>  +  Aa)sin0+h  (8) 

* 

y  =  a  cos  (  cj>  4*  A<f> )  sin0  4-  b  sin  (  <J>  +  Aa  )  cos  0  4-  k  .  (9) 

Using  the  double  angle  formulas, 

cos  (  $  4-Aa  )  =  cos^cosA a  -  sin^sinAa  (10) 

sin  (  <(>  +  A cj>  )  -  coscj)SinA(J)  +  sincjjcos  a cj>  ,  (11) 


we  need  only  calculate  the  cosines  and  sines  once  for  each  ellipse  drawn. 
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A  plotting  program  can  be  written  as  follows: 

Calculate : 

a  sin  0  ,  a  cos  0  ,  b  sin  0  ,  b  cos  0  ,  sinAtj),  cos  A<j) 

Initialize: 

Pen  up 

cos  <|>  »  1  and  sincj)  =  0 

Move  pen  to  (a  cos  0  4-  h,  a  sin  0  4*  k) 

Pen  down 
Then  do: 

Using  the  present  values  of  cos  <j>  and  sin  (j>  ,  calculate  the  next 
cos  b  and  sin  4>  using  Equations  (10)  and  (11),  then  update  cos 
4>  and  sin  $  ,  and  then  draw  to  the  new  location  defined  by 
Equations  (8)  and  (9). 

Loop  until  done. 

Within  the  Do  loop  there  are  only  8  multiplications,  4  adds,  and  2  subtractions. 
Note  that  no  trigonometric  function  need  to  be  calculated  in  the  loop. 

Figure  3  is  a  sample  program  using  the  above  algorithm  and  Figure  4  is  a  plot  of 
the  resultant  data. 


11 
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PROGRAM  ELLIPSE 

THIS  PROGRAM  IS  AN  EXAMPLE  OF  PLOTTING  AN  ELLIPSE  USING  THE 
ALGORITHM  FROM  CHAPTER  3  OF  THIS  REPORT. 

THE  ELLIPSE  HAS  AN  INCLINATION  OF  45  DEGREES,  A  SEMI -MAJOR 
AXIS  OF  1,  A  SEMI -MINOR  AXIS  OF  0.5  AND  IS  CENTERED  AT  (1,0.5). 
THE  DELTA  PHI  INCREMENT  IS  1.8  DEGREES. 

INITIALIZE  THE  GIVEN  PARAMETERS . 

PI  -  3.14159 

DELTA  -  45 . 0*PI/180 . 

DPHI  -  1 . 8*PI/180 

A  =  1.0 

B  =0.5 
XH  =  1.0 

YK  =0.5 

CALCULATE  THE  SINES  AND  COSINES  REQUIRED. 

SD  =  SIN(DELTA) 

CD  =  COS (DELTA) 

SDP  =  SIN(DPHI) 

CDP  =  COS (DPHI) 

CALCULATE  A*SIN(DELTA) ,  A*COS (DELTA) ,  B*S IN (DELTA)  AND  B*COS (DELTA) . 
ASD  =  A*SD 
ACD  -  A*CD 
BSD  =  B*SD 
BCD  =  B*CD 

INITIALIZE  THE  REQUIRED  VARIABLES. 

CP  =1.0 

SP  =0.0 

THE  PEN  IS  UP  AT  THIS  POINT;  I.E.,  NOTHING  HAS  BEEN  WRITTEN  YET. 

MOVE  THE  PEN. 

X  =  ACD+XH 
Y  -  ASD+YK 

PUT  THE  PEN  DOWN;  I.E.,  WRITE  THE  DATA  POINTS  FOR  LATER  PLOTTING. 
WRITE(10,*)  X, Y 

THIS  IS  THE  BEGINNING  OF  THE  LOOP.  WE  WILL  LOOP  UNTIL  DELTA  PHI 
COVERS  360  DEGEES;  I.E.,  1.8*200. 

DO  100  1=1,200 
C  CALCULATE  EQUATIONS  10  AND  11. 

CPN  =  CP*CDP  -  SP*SDP 
SPN  =  CP*SDP  +  SP*CDP 
C  UPDATE  COS (PHI)  AND  SIN (PHI) . 

CP  =  CPN 
SP  =  SPN 

C  DRAW  TO  THE  NEW  LOCATION  DEFINED  BY  EQUATIONS  8  AND  9. 

X  =  CP*ACD  -  SP*BSD  +  XH 
Y  =  CP*ASD  +  SP*BCD  +  YK 

C  THIS  WRITE  STATEMENT  IS  OUR  RECORD  OF  THE  DRAW. 

WRITE (10 ,*)  X,Y 
100  CONTINUE 
STOP 
END 


FIGURE  3.  AN  EXAMPLE  PROGRAM  FOR  PLOTTING  AN  ELLIPSE 
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PLOT  OF  THE  ELLIPSE 


♦ 


FIGURE  4 


A  PLOT  OF  THE  DATA  FROM  THE  ELLIPSE  PROGRAM 
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